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On the statistical mechanics of prion diseases 
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(February 1, 2008) 

We simulate a two-dimensional, lattice based, protein-level statistical mechanical model for prion 
diseases (e.g.. Mad Cow disease) with concommitant prion protein misfolding and aggregation. Our 
simulations lead us to the hypothesis that the observed broad incubation time distribution in epi- 
demiological data reflect fluctuation dominated growth seeded by a few nanometer scale aggregates, 
while much narrower incubation time distributions for innoculated lab animals arise from statistical 
self averaging. We model 'species barriers' to prion infection and assess a related treatment protocol. 



PACS Indices: 82.30.-k, 82.35.Pq, 87.15. Cc 

Transmissible neurodegenerative prion diseases, such 
as Mad-Cow Disease (BSE) and Creutzfeldt-Jakob Dis- 
ease (CJD) in humans, increasingly represent a serious 
pubhc health threat 0] . Prusiner and collaborators have 
shown that the infectious agent (prion) in these diseases 
consists of a quantity of a misfolded form (PrP'^'^) of the 
^200 amino acid PrP*^ protein which is expressed ubiq- 
uitously in mammalian neurons j^]. The PrP*^ proteins 
normally reside on neuron surfaces , and the more hy- 
drophobic PrP'S'c forms tend to aggregate: large (micron 
scale) PrP'^'^ amyloid plaques are a common post-mortem 
feature of brain tissues, and fibrils are observed in in vitro 
experiments. Nucleic acid free propagation demands that 
Pj.pSc autocatalyze their own formation by helping to 
convert more PrP'^ into PrP-^'^ [HH]- Given that mam- 
malian PrP'^ differ by only 5-10% in amino acid composi- 
tion and that variant CJD correlated with BSE has been 
observed in England, the efficacy of "species barriers" in 
limiting transmission is of considerable interest. 

Several facts suggest that prion disease may be driven 
less by complex biology (as in Alzheimers' disease) than 
by physico-chemical processes, including: (1) The "uni- 
versal character" of sporadic CJD, observed globally 
without evident spatio-temporal clustering at an annual 
background death rate of 0.5-1.5/million and mean onset 
age of about 63 years Q. (2) Highly reproducible incu- 
bation time vs. infectious dose distributions observed in 
laboratory animal studies, with a power law relation be- 
tween mean incubation time and dose concentration D 
0. (3) The close temporal proximity of disease onset 
and death. (4) PrP-^c may be grown in vitro without 
biological processing and is toxic to neurons, though not 
yet demonstrably infectious |^ . The following question is 
thus raised: is the incubation time of prion disease dom- 
inated by a fundamental physico-chemical nucleation and 
growth of PrP^'^ protein aggregates ? 

We present simulation results for a simple two- 
dimensional, lattice based statistical mechanical model 
for prion disease based at the protein level. It is all but 
hopeless to develop an atomic level simulation capable of 
spanning the 21 temporal orders of magnitude between 
picosecond scale intra-protein motion and the potentially 



decades long incubation times of large prion aggregates. 
Our simulation is a bridge between the short distance 
and time scales covered in individual protein models and 
the long time, macroscopic realm of chemical kinetics. 
We use our simulations to identify the protectorate of 
principles necessary to describe the aforementioned uni- 
versal features of prion disease. Specifically, we find: (1) 
Concommitant PrP'^'^ autocatalysis and stable aggrega- 
tion can explain the long disease incubation times, and 
the difference between sporadic (unseeded) and infectious 
(seeded) onset. This principle is consistent with data for 
yeast prions [|l^, is closely related to the earlier nucle- 
ation theory of Lansbury , and avoids the need to fine 
tune model parameters with separated autocatalysis and 
aggregation (2) The distributions of disease incu- 
bation times (well defined onset edge, broad fluctuation 
driven spectrum for dilute prion doses, narrow short time 
shifted distribution for more concentrated prion doses) 
stongly suggest a central role for growth from a minimally 
nanometer scale seed (with order 10 PrP'^'^ proteins). 
With these assumptions we obtain a one parameter fit 
to the BSE incubation time data, an order of magnitude 
correct estimate of the mean incubation time for BSE and 
infectious CJD, and we argue that self-averaging explains 
a narrowing of incubation time distributions observed in 
the laboratory [@jll) . Finally, by considering both stick- 
ing and interconversion rate differences for prions of dif- 
ferent species, we provide a theoretical underpinning to 
data for the species barrier efhcacy and to proposed treat- 
ments based upon coating of incipient prion plaques. 

Our assumption of two-dimensionality is consistent 
with the observation that PrP"^ resides on neuron surfaces 
and autocatalysis and aggregation must certainly initiate 
there. The lattice structure should be irrelevant at long 
times. By assuming a triangular lattice we minimize the 
anisotropy (which appears first in rank six tensors for 
hexagonal lattices). This simulation yields the crucial 
qualitative features necessary to specify the "prion pro- 
tectorate." Our resulting prion aggregates are essentially 
amorphous, rather than quasi-one-dimensional filaments 
as found in some in vitro experiments. We shall discuss 
an anisotropic simulation yielding fibrils elsewhere. 



Our model is based upon the energy landscape schema- 
tized in Fig. 1. Individual lattice cells take on three val- 
ues: unoccupied (implicitly filled with water), occupied 
by a PrP"^ protein, or occupied by a PrP'^'^ protein. At 
the beginning of each simulation, we choose a monomer 
PrP"^ configuration distributed randomly, with a number 
(0-4) of PrP'S'c seeds located randomly as well. 

In the spirit of cellular and molecular automata, we es- 
chew an explicit Hamiltonian based molecular dynamics 
in favor of the following updating rules, specified in time 
units of one full lattice sweep (d is the nearest neighbor 
coordination number): 

(1) Identify each occupied cell as PrP^ or PrP^'^ by 
checking its number of nearest neighbors. 

(2) For each Prl^ , move one step in a random di- 
rection if that adjoining site is unoccupied. 

(3) Identify aggregates (size Np) by PrP^'^ presence; 
move them per Rule (2) with probability 1 / yjNp . 

(4) For d = 0, 1 the protein remains PrP^ . 

(5) For d = 2, the protein is PrP^ or PrP^'^ with equal 
probability 

(6) For d > 3 the protein is PrP^'^. 

Rules (1-3) ensure diffusive motion of PrP*^ monomers 
and PrP'^'^ aggregates. Rules (4-6) specify protein in- 
teractions and autocatalysis, and reflect the relative hy- 
drophobicity of PrP'^'^. To assure the biologically plausi- 
ble condition of constant average PrP*-^ concentration, at 
the start of a new sweep we randomly place a new PrP*^ 
for each lost in the previous update. 

This procedure is repeated for up to millions of full 
sweeps through the lattice. Note that a sweep defines 
the minimum time scale for conversion of PrP*-^ to PrP'^'^, 
which is presumably of order a second or less in real time. 
We work with maximum lattice sizes of = 4 x 10''. 
Practically, we are able to run at areal monomer con- 
centrations c > 0.1%, which are likely to be about three 
orders of magnitude larger than in the brain. 

We run our seeded simulations repeatedly following 
the largest nano-aggregates until they reach a size of up 
to 0.4% of the lattice. Each seed consists of 10 PrP'^'^ 
monomers [T^ j For a given aggregate size, we have stud- 
ied the properties of the distribution of lifetimes to reach 
that size. For a final aggregate size of ^=0.2%7V (80 
monomers) we display the corresponding distributions 
Pd in Fig. 2(a) for several values of c and D = 1. Pd 
displays greater fluctuations for small c, and increasing 
c shifts Pd to short times. Fluctuations dominate the 
growth from the small seed aggregate. For larger A we 
find that the onset time grows sublinearly in A, while 
the distribution width grows weakly with A. To within 
fluctuations, the Pd curves for different c collapse when 
scaled by the mean time tm- (The curves broaden some- 
what for D ~ 2, 3, 4.) This implies the scaling law 
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In the range of concentrations studied, we can fit t,m{c) 
well by either: (a) the polynomial form im(c) ~ Ac^^ + 
Bc^^ (the first term refiecting monomer aggregation, the 
second representing dimer aggregation), or (b) a power 
law form, t„i{c) ~ c-"(^), with a{D) = 1.66(£) = 
1), 1.49(1? = 2). For finite A and c — > the exponent 
a{D) should tend to two, refiecting dominant dimer ag- 
gregation. Fit (b) is consistent with a nontrivial scaling 
limit as .4 — + oo, c ^ 0. Both fits to t,n{c) give order 
of magnitude equivalent extrapolations to lower c rele- 
vant to biological contexts. Assuming im(c) ^ at 
fixed c, we estimate f3 < I with large numerical uncer- 
tainty. Taken together, these data lead us to a conjec- 
ture of asymptotic compression of Pd'- for c ^ 0, an in- 
creasingly large range of A will have essentially the same 
scaled shape, because the time t2 it takes to go from A 
to say, 2A will be much smaller than tm for going from 
10 ^ A. 

We believe that other processes, such as protease at- 
tack or competition between PrpS'c.pj.pSc ^^^^ pj.pSc_ 

neuron binding will limit prion aggregation on the neu- 
ral surface and cause fissioning and spreading of nano- 
aggregates to other neurons. Such fission is necessary for 
exponential growth of infectious agents, as noted previ- 
ously [Q. At the order-of- magnitude level t2 sets the 
doubling time for disease spread, assuming that inter- 
cellular diffusion rates exceed intra-cellular ones. Given 
t2 « tm (at biological c values) both in our fission hy- 
pothesis and experiment, slow growth fluctuations for 
small seed aggregates can dominate the overall incuba- 
tion time for disease. For example, in hamster experi- 
ments the first observable levels of PrP'^'^ in brain tissues 
occur at about 90 days, with symptom onset at 120 days, 
and a doubling time of 2 days [|l^. Given the conjec- 
tured asymptotic compression, the incubation time dis- 
tribution will be relatively insensitive to a wide range of 
nano- aggregate sizes for fission parents or progeny. 

We stress several points of agreement with observa- 
tion and laboratory data. First, we always find death 
(runaway growth of a killer aggregate) in the presence 
of infection, in agreement with experiment and the ag- 
gregate nucleation hypothesis for prion disease. Second, 
our rough scaling of tm with D agrees qualitatively with 
the extraordinarily reproducible dose vs. incubation time 
data of Prusiner and collaborators for prion infection in 
laboratory mice, for which r/ ~ Third, our life- 

times are in order of magnitude agreement with observed 
BSE incubation times and new variant CJD death ages, 
which we can prove by exploiting our concentration scal- 
ing. For a single seed, at our lowest concentration of 0.1% 
our distribution peaks at about 10^ sweeps. Assuming a 
literature estimate of concentration of in vivo prions in 
solution of 100 nanomole/liter [Esl and a protein diame- 



ter of order one nanometer suggests a dimensionless areal 
concentration of about 10~'^%. With the above concen- 
tration scahng this gives an incubation time of about 10^ 
sweeps. If we take the 5 year mean incubation time for 
BSE this imphes our basic misfolding time (one sweep) 
is about one second, which is reasonable at the order 
of magnitude leveL Finally, by simply scaling time for 

= 80 by the mean time tm, we achieve an exceptionally 
good fit (A) of the inferred incubation time distribution 
for BSE in the United Kingdom (for cattle born in 1987) 
p6[ to our scaling curve as shown in Fig. 2(b). This 
supports our hypothesis that runaway autocatalytic ag- 
gregation seeded by a few prion nano-aggregates triggers 
the disease. This is very plausible since (i) infection is 
believed to have come from rendering plant offal derived 
from many animals after considerable processing, leading 
to substantial dilution, and (ii) prion transmission by 
oral consumption is estimated to be less efficient than 
innoculation by a factor of 10^. 

Sporadic CJD develops in the absence of any infec- 
tion or genetic predisposition to prion disease. We have 
studied the lifetime distributions in this case and find 
that they are very different in shape than the seeded 
runs, being much broader and dimer fluctuation domi- 
nated (peaks scale as c~^). For our estimated biological 
concentration of c = 10~^%, this could give a mean onset 
time as large as 10'^ years. Assuming a constant height 
tail from the minimum single seed onset to the mean hu- 
man lifespan taken with the 1 part in 10^ odds of dying 
from sporadic CJD gives an estimate of a total death 
probability of a part in 10^ up to this mean on-set time. 
This is consistent with our simulations, but not resolv- 
able currently. 

In contrast to the breadth of our single seed lifetime 
distributions, in the event of dosing by a greater number 
of seeds, self-averaging will narrow the distribution, col- 
lapsing it towards the onset (shorter times). Assuming 
growth proceeds independently for different seeds (which 
is quite reasonable for the low concentrations in living or- 
ganisms) distributions would narrow by y/D. Several in 
vivo experiments show distribution time narrowing with 
increased doseage, but are not amenable to quantitative 
analysis at this time 

We turn now to the species barrier. To simulate inter- 
species transmission we replace the integer-valued coor- 
dination, the key parameter which determined the stabil- 
ity of different protein conformations in the earlier stud- 
ies, by a new variable x ~ N + N' x P for PrP'-^ and 
x' ^ N' + N X P', for PrP*^', where N is the number of 
Pj.pC,Sc neighbors and N' the number of PrpC.-S'c' neigh- 
bors. The parameters P and P' (chosen less than unity) 
provide a measure of the reduced effectiveness of inter- 
species conversion. To deal with the continuous variable 
X, rather than the integer coordination used earlier, we 
represent the PrP'-^ to PrP'^'^ conversion probability f{x) 
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with a similar definition for PrP*^' -^PrP"^'^' conversion. 
A separate parameter describes physisorption to an ag- 
gregate of alien prions. We allow a PrP^^'Sc (jiyg^^gj- Qf 
mass M (M — 1 being a monomer) to desorb from an 
alien prion aggregate with probability /\/M, where 
q is the physisorption parameter, and K is the number 
of units in the cluster adjacent to the alien aggregate. 
We define q' for pj-pC.'S'c/ gjniilarly. Thus, we describe 
species barriers in terms of the two parameters affecting 
conversion rates (P, P') and physisorption {q,q'). Mak- 
ing P ^ P' and q ^ q' allows for an asymmetry in infec- 
tivity, and we can model one of the most striking aspects 
of the species barrier, namely its asymmetry {e.g. mice 
infect hamsters well, hamsters infect mice poorly). 

We have carried out several simulations relevant to the 
species barrier. First, we have placed alien prion 'walls' 
on a boundary, representing a single large aggregate from 
another species. With suitable choice of parameters, 
this leads to the asymmetric species barrier between two 
species. Namely, disease onset in one species (with large 
aggregates) and only sub-critical PrP-^'^ concentration 
near the wall in the other. Recent experimental data 
has also shown the build up of such sub-critical PrP'^'^ 
concentration in the asymmetric interspecies transmis- 
sion III]. 

Second, suppose the asymmetry is such that P << P'. 
In this case, PrpS'c jg niore favorably formed from PrP"-^' 
in the presence of a seed aggregate of PrP'^'^ than PrP'^'^ 
is formed from PrP'-^ in the presence of a PrP^'c geed. 
Thus, by injecting a large enough initial dose of alien 
Pj.pC/ proteins into the organism, it should be possible 
for PrP*^' to compete favorably with PrP*-^ for aggrega- 
tion, and thus extend the incubation time. Such proto- 
cols have been tested experimentally: (i) Experiments 
with the coating dye molecule Congo Red reveal 
a surprising non-monotonic dose vs. lifetime relation: 
small Congo Red concentrations yield a reduced time for 
incubation, while larger Congo red concentrations signif- 
icantly extend the incubation time, (ii) In vitro exper- 
iments show that the above coating scenario with alien 
prions works p9[ | when the initial PrP*^' concentration 
c' > c. 

We find both of these experimental features in our sim- 
ulation as shown in Fig. 3. For c'/c < 2, the incuba- 
tion time is slightly shortened, while for c'/c > 2 the 
incubation time is increased. We find that the lifetime 
shortening arises because once a few PrP'-^' aggregate, 
they partially block the motion of adjacent PrP'-^ pro- 
teins, enhancing the misfolding probability of the latter. 
However, further PrP*^ misfolding is significantly reduced 
once PrP'^'^' coats the seed. 

A crucial test of our theory would be to monitor, in 
vitro, the lifetime distributions of normal prion proteins 



seeded by small aggregates as a function of both c and 
D to look for the scaling behavior. If it is not possible 
to easily achieve D = 1,2, presumably the distribution 
of Fig. 2(b) could be estimated by extrapolation from 
successively diluted PrP'^'^ doses. This would also test 
the extent to which the self-averaging applies. If our 
single seed hypothesis is correct, it is important to assess 
how passage occurs upon infection. 

Theoretically, it is important to extend our model by 
including internal degrees of freedom for the proteins 
which can provide an orientational dependence to the ag- 
gregation, allowing for linear fibril formation as opposed 
to amorphous aggregates. This connects naturally with 
simulations on individual proteins, and also allows for 
the study of different protectorates for collective protein 
phenomena. 
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FIG. 1. Schematic energy landscapes for model prion pro- 
teins. Upper figures show proteins in varying coordination 
environments corresponding to rules (1-6), lower figures show 
schematic energy landscapes. 
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FIG. 2. Aggregation lifetime distributions PD{t,c) for a 
single seed (D—l). Seeds have ten proteins, simulations ter- 
minate at aggregate size ^ = 80, time is measured in units 
of 10"* full lattice sweeps, (a) Dependence upon areal con- 
centration c of normal prion proteins, (b) Scaled distribu- 
tion Fi{t/tm) vs. t/tm- tm is the mean lifetime for a given 
concentration. Points are inferred incubation time data for 
BSE-infected cattle in the United Kingdom born in 1987. 
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FIG. 3. Mean incubation time tm vs. relative concentra- 
tion of alien prions cf/c (test of coating treatment; see text 
for discussion). 



